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A phenomenological Lagrangian approach is employed to study the electromagnetic properties of 
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construct a general form for the electromagnetic one- and two-body transition operators formulated 
in terms of the nucleon fields, which are then used in the calculation of the electromagnetic deuteron 
form factors. One of the two-body operators is responsible for explaining the quadrupole moment 
form factor. We show that in our approach the data on the deuteron form factors as well as on the 
differential cross section of elastic electron-deuteron scattering are well explained. 
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I. INTRODUCTION 



The study of the electromagnetic properties of the deuteron has a long and rich history (for some recent reviews, 
see, e.g. The deuteron, as a spin-1 particle, is usually beheved to be a weakly bound system of a proton and 

a neutron with a binding energy ejj ~ 2.22 MeV. Since the electromagnetic (EM) properties of the deuteron can also 
shed light on the EM form factors of the neutron as well as on nuclear effects on the form factors, the study of the 
deuteron form factors with a EM probe is of great interest. 

The matrix element for elastic electron-deuteron (eD) scattering in the one-photon approximation is 



M = -^Ue{k'h^Ue{k)J^ip,p') (1) 



where k and k' are the four-momenta of initial and final electrons and J7^(p,p') is the deuteron EM current - 

J^'^iP,?') - - (Gi(g2)e'* • e - ' qe'* • 9) (p + p')^ - G^iQ^) {e^e'* ■ q - e'^e ■ , (2) 

where mu is the deuteron mass, e(e') and p{p') are polarization and four-momentum of the initial (final) deuteron 
with q — p' — p being the momentum transfer. The three EM form factors Gi,2,3 of the deuteron are related to the 
charge Gc, quadrupole Gq and magnetic Gm form factors by 

2 

Gc=Gi + -tGq, Gm — G2 , Gq = Gi - G2 + {I + t)G3, t — . (3) 



These form factors are normalized at zero recoil as 



Gc(0) = l, GQ{0)^mlQD^ 25.83, Ga/(0) = —^ic = 1.714 , (4) 

where rriN is the nucleon mass, Qd and fj.D are the quadrupole and magnetic moments of the deuteron. Since the 
deuteron is a spin-1 particle it has three EM form factors in the one-photon-exchange (OPE) approximation, due 
to current conservation and the P and C invariance of the EM interaction. The three form factors Ge.m,q can be 
determined by measuring the unpolarized, elastic eD differential cross sections and one of polarization observables, 
like the deuteron polarization tensor 

T20 = {^ItGcGq + Ir'Gl +^-^1 + 2(1 + r)tan2|) G^) (5) 

with S = A{Q^) + B{Q'^) tan2(6'e/2). The two form factors A{Q^) and B{Q^) are related to the EM form factors of 
the deuteron as 

A^GI + ^tGI + ^t'GI, B^^t{1 + t)GI,. (6a) 

According to the Rosenbluth separation Q , the elastic scattering of an unpolarized electron from the deuteron results 
in an 0{a^) differential cross section da/d£l = {da / dVl)yiott S, where 9^ is the electron scattering angle in the 
laboratory frame of the collision and {da/dfl)Mott is the Mott cross section. 

The theoretical study of eD elastic scattering and deuteron EM form factors has been performed in different ap- 
proaches potential models, phenomenological models including quark, meson and nucleon degrees of freedom, 
effective field theories, etc. (for an overview see [l|-[3)- In the present work we apply a phenomenological Lagrangian 
approach to study the EM form factors of the deuteron. We consider the deuteron as a weakly bound system of 
proton and neutron. The coupling of the deuteron to its constituents is determined by the compositeness condition 
Z = [IB, [13] I which implies that the renormalization constant of the hadron wave function is set equal to zero. 
Note, that this condition was originally also applied to the study of the deuteron as a bound state of proton and 
neutron [l6j . Then it was extensively used in low-energy hadron phenomenology as the master equation for the 
treatment of mesons and baryons as bound states of light and heavy constituent quarks (see Refs. [11]). In Refs. [l^ 
this condition was used in the application to hadronic molecule configurations of light and heavy mesons. 

To study the EM form factors of the deuteron, in a first step we employ the empirical EM couplings of the photon to 
the proton (or neutron) in the one-body operators to set up the well determined part of the photon-deuteron coupling. 
As a consequence of the non-local description of the deuteron bound state direct contact terms are then included to 
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guarantee local gauge invariance of the EM interaction. In the last step we introduce additional, phenomenological 
two-body operators, which are assumed to represent e.g. meson-exchange currents and in turn imply a D-wave 
component in the deuteron wave function. Parameters of two-body operators are deduced from a fit to the available 
data. 

The paper is organized as follows. In Sec. II, we discuss the basic notions of our approach: the coupling of the 
deuteron to its constituents involving the compositeness condition and the derivation of the EM one- and two-body 
operators contributing to the form factors of the deuteron. In Sec. Ill we discuss the numerical results for the deuteron 
EM form factors as well as the two form factors A(Q^) and B{Q^) entering in the differential cross section of the eD 
elastic scattering. In Sec. IV we give our conclusions. 



II. APPROACH 



A. Deuteron as a proton— neutron bound state 



In this section we discuss the formalism for the study of the deuteron interpreted as a hadronic molecule - a weakly 
bound state of proton and neutron: \D) = \pn). We write the deuteron mass in the form 1710 = 2mN — en, where 
niN = mp = 0.93827 GeV is the nucleon mass and cd — 2.22 MeV is the binding energy. Based on our approach, the 
coupling of the deuteron to its two constituents - proton and neutron, is 

Cd{x) = cjdDUx) f dy^D{y^)p{x + v/2)Crn{x - v/2) + H.c, (7) 



where C = 7^7^ is the charge conjugation matrix and x is the center-of-mass (CM) coordinate. The correlation 
function $d characterizes the finite size of the deuteron as a pn bound state and depends on the relative Jacobi 
coordinate y. A basic requirement for the choice of an explicit form of the correlation function is that its Fourier 
transform vanishes sufRciently fast in the ultraviolet region of Euclidean space to render the Fcynman diagrams 
ultraviolet finite. We adopt the Gaussian form, ^d{p\,) = cxp(— p|j/A|,) , for the Fourier transform of the vertex 
function, where is the Euclidean Jacobi momentum. Here, is a size parameter that characterizes the distribution 
of the constituents in the deuteron. 

The coupling constant (/d in Eq. ([7]) is determined by the compositeness condition [l6l.[l7j. which implies that the 
renormalization constant of the hadron wave function is set equal to zero: 

Zd = l-S'^(mi) =0. (8) 

Here, I]'^(m|,) = g^^'D{''^ii) the derivative of the transverse part of the mass operator , conventionally split 
into the transverse E/j and longitudinal pa-rts as: 

^'^f{p)^gfj:n{p') + ^^Up'), (9) 

where 5"^ = g"^ — p"p^ /p"^ , g'^Pa = . The mass operator of the deuteron is described in Fig.l. For a fixed value 
of the size parameter A^i the coupling is determined according to the compositeness condition. 

1/2 

To clarify the physical meaning of the compositeness condition, we reiterate that the renormalization constant Z^j 
can also be interpreted as the matrix element between the physical and the corresponding bare states. For Zjj = 
it then follows that the physical state does not contain the bare one and hence the deuteron is described as a bound 
state of the proton and neutron. As a result of the interaction of the deuteron with its constituents, the deuteron is 
dressed, i.e. its mass and its wave function have to be renormalized. 

Following Eq. ^ the coupling constant go can be expressed in the form: 

1 ^ ^00 .00 ^^^^ / a + P + AaP 



r dad(3 ( o( j^R\ ^ a + (3 + Aa(3 ^\ 

L (l + a + /3)3 ^-<-^(" + ^)^^+ 2(l + . + ;3) H 



(l + 2a)(l + 2/3) 2 \ (l + 2a)(l + 2/3 ) 
2(l + a + /3) ■ 4(l + a + /3)2 ''^j ^ 2(1 + a + /3) 



X I (« + /? + AaP) + , ^ , , ^ , o\2 ^''D] + 0/1 , -. , ^^ 1' (10) 



where = mH / with H ~ N^D. 
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B. Matrix element of the photon-deuteron interaction 



To calculate the deuteron EM form factors, we construct the electromagnetic transition operator including one- and 
two-body parts and formulated in terms of nucleon degrees of freedom - the constituents of the deuteron. Note, that 
the direct coupling of the deuteron with the photon field vanishes because the Zj^ factor equals zero. It guarantees 
that double counting is avoided. 

We construct the one- and two-body operators in a general phcnomenological form, which principally includes all 
possible corrections (e.g. meson-cloud effects) and the dependence on the photon momentum (form factors). The 

one-body operator (or its Fourier transform) with J^^\q) = Jfi^^ig) + Jj^^^ {o) contains two terms. The first term, 
J^^^ (q), is generated by the coupling of nucleons to the electromagnetic field: 



j;^^(g) = ] d^:,e~^^^N{x)[^rF-{<i') + ::^Fi\q')^N{x), (11) 

where F-^ (q^) and F^{q'^) with N — p,n are the conventional Dirac and Pauli form factors of the nucleon [see 
Fig. 2(a)]. The second term, J^^^{q), is generated by gauging the nonlocal strong Lagrangian Co [see Fig. 2(b)]. To 
restore electromagnetic gauge invariance in Co, the proton field should be multiplied by the gauge field exponential 
(see further details in p^): 

X 

p{y) ^ e-^(-^^^^)p(y) , I{x, y,P)= f dz,A''{z). (12) 



For the derivative of I{x,y,P) we use the path-independent prescription suggested in [201] which in turn states 
that the derivative of /(x, y, P) does not depend on the path P originally used in the definition. The nonminimal 
substitution (|12p is therefore completely equivalent to the minimal prescription. Expanding the exponential e**^^*^^'^'^' 
in powers of the electromagnetic field and keeping the linear term (corresponding to the vertex D''pn^ + H.c.) we 
generate an additional contribution [see Fig. 2(b)] to the electromagnetic one-body operator: 

x+y/2 

jDNN^^^ = ^tgo I d^xd^yDl{x)^D{y^)p{x + yl2)CYn{x - y/2) f dz^e'^'^' + H.c. (13) 









= J d'^xe' 


4^-^\q) 


~ j d'^xe" 




— J d^xe~ 



X 

There is a number of possible contributions to the two-body operator jj^\q) — J^'^{q). We restrict to the three 
simplest terms with the smallest number of derivatives [see the general diagram of Fig. 2(c)]: 

(14a) 

'''^^giF,^''{q^)n{x)rCp{x)p{x)Cj^iaf,,q''n{x) + H.c. , (14b) 

-^^-g^Fi"'{q^)n{x)fiCp{x)p{x)Cta^,q''n{x) + H.c. , (14c) 

^^^^{q')[n{x)rCp{xMd,, - d^Mx)C^Mx)] , (14d) 

where gi and F,[^^{q'^) are the phcnomenological electromagnetic two-body nucleon couplings and form factors, 
respectively. 

To calculate the EM form factors of the deuteron we project the dressed operator Jfj,{q) = J^p (q) + j/i'^' {q) between 
the deuteron states: 

{D{p')\J,{q)\D{p)) = (27r)4 5\p' - p - q) J^{p.p') , (15) 

where J^(j)^p'^ is the deuteron EM current given by the expression ([2]). The diagrams contributing to the matrix 
element (fTSj) are shown in Fig. 3 - the diagrams generated by the one-body currents [Figs.3(a)-3(c)] and by the two- 
body currents [Fig. 3(d)]. To evaluate these diagrams we take the T-product of the EM current in the ^-matrix 
defined in terms of the interaction Lagrangian Ld and use the standard free fermion propagators for the nucleons in 
the loops. Let us stress that a similar approach to the one presented here was previously developed in [T5| . However, 
in Rcf. 15] the authors did not consider two-body operators. Also note that the deuteron EM current generated by 
both the one-body and two-body nucleon operators is manifestly Lorentz covariant and gauge invariant. 
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We want to point out again that our approach is purely phenomenological in the sense that we do not calculate the 
one- and two-body operators from microscopic models, but constrain their forms using a justified physical background: 
in particular, the two-body operators (fT4| introduced in our considerations have an explicit nonrelativistic limit. They 
correspond to the ones generated in the context of the effective field theory description [3]. 

The correlation function $£i of Eq. ([7]), modeling the distribution of nucleons in the deuteron, has in the full formal- 
ism no direct connection to the quantum-mechanical wave function of the deuteron. Comparison of our formalism to 
the nonrelativistic approaches, dealing with the deuteron wave function, can be performed only on the level of matrix 
elements. In this vein the expectation value of the two-body operator Jff^'^{q) between deuteron states incorporates 
the D-wave admixture to the deuteron wave function in the context of potential models. This admixture is necessary 
to explain the quadrupole moment of the deuteron. In the context of effective field theory [l3| it was shown before 
that the nonrelativistic analog of the Jff^''^{q) operator explains the quadrupole moment/form factor of the deuteron. 
Therefore, we introduce this operator in consistency with previous observations. 



III. NUMERICAL RESULTS 



To describe the EM form factors of the deuteron we have the following input: the size parameter Ad, describing 
the distribution of nucleons in the deuteron, the EM form factors of the nucleons Fi'2{<f') and a set of parameters in 
the two-body operators - couplings constants gi and form factors F^-^ {q^). In fixing the value of the parameter A^ 
we use the following constraint. In the non-relativistic approximation the vertex function ^oi—p^) represents the 
wave function of the deuteron. Therefore, a constraint condition for A^ is set by the deuteron size that, according 
to potential model calculations, is bound as d?""^!) < 0.02 GeV^ Employing this condition, we expect that A^ is 
less than 0.5 GeV. 

For the set of the EM nucleon form factors -Ff';^(9^) we use parametrizations in a fit to the corresponding data. 
We use two available forms - the Mergell-Meissner-Drechsel (MMD) parametrization [2H and the Kelly parametriza- 
tion [23. We stress that when we restrict to the use of one-body transition operators only we cannot reproduce the 
data on the EM form factors of the deuteron. In particular, we cannot correctly predict the Q^-dependence of the 
deuteron form factors - the nodes of the charge and magnetic form factors at alsout 0.7 GeV^ and 2 GeV^ [2^ are 
not reproduced. Also, the quadrupole form factor cannot be explained - with the one-body operators present only 
its normalization is completely underestimated. As is known from potential models the D-wave component in 
the deuteron wave function mainly contributes to the quadrupole moment and the quadrupole moment vanishes if 
only the S-wave deuteron wave function is considered. Because at this level our correlation function P^) does 

not contain an explicit admixture of a I?- wave component, the result for the quadrupole moment with one-body 
operators is negligibly small when compared to data. 

In the present approach we propose and test if a suitable choice for the additional two-body operators can reproduce 
the full quantitative structure of the deuteron form factors. In particular, as will be shown, the two-body operator 
Jfi^'^iq) will give the dominant contribution on top of the one-body structures to fully explain the quadrupole 
form factor. We proceed by fixing the related parameters g2 and F^^{q'^) of Jft^''^{q) to reproduce the measured 
quadrupole form factor of the deuteron. The couplings 31,3 and the form factors F^^ (Q^) of the additional two-body 
operators are adjusted to obtain a refined fit to the charge and magnetic form factors, including the nodes at about 
0.7 and 2 GeV^ [23|. For the form factors F^^^{Q'^) we use the following parametrization: 

where A, Ai and A2 are the size parameters (for Fi^iCf') and F^^ [Q"^) we use the same ones). The form factors 
Fi^iCp) and F^^ [Q^) should vanish at zero recoil such that the deuteron charge is not renormalized (to preserve 
the charge conservation), while the form factor F^^{Q'^) does not affect current conservation and, therefore, does 
not vanish at zero recoil. In Table I we present the results for the fit parameters 51,2,3, A and Ai,2 using the data on 
the deuteron EM form factors for two different parametrizations [21I. [2^ of the nucleon Dirac and Pauli form factors. 

Our numerical results for the deuteron charge, magnetic and quadrupole form factors as well as for the two form 
factors, A{Q'^) and B{Q^), entering in the differential cross section and for the polarization tensor 120 are shown in 
Figs. 4-18. In particular, in Figs. 4-6 we present the results of our approach for the Gc{Q^), Gq{Q^) and Gm{Q^) 
form factors (total contribution of the set of the diagrams in Fig. 3). In Figs. 7-9 we analyze the separate contributions 
of the graphs of Fig. 3 to G'c(Q^), Gq(Q^) and Gm{Q^)- In Figs. 10-15 we show the separate contributions of the 
two-body operators generating the diagram of Fig. 3d on a logarithmic scale for |Gc((5^)|, \Gq{Q'^)\ and \Gm{Q'^)\ 
(Figs.lO, 12, 14) and on a linear scale for Gc{Q^), Gq{Q^) and Gm{Q^) (Figs.ll, 13, 15). In Figs.l6 and 17 we 
plot our results for the form factors A{Q^) and B{Q^). Finally, in Fig. 18 we present the results for the deuteron 
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polarization tensor T2q{Q^). We select 9e = 70° to make a comparison consistent with the experimental data. We also 
perform a comparison with other theoretical calculations. We present our results for the two sets of parametrizations 
of the nucleon EM form factors - the MMD [2l| and the Kelly parametrization f23| . We also compare our fit to the 
parametrization derived in Ref. ^2§\ (denoted as TGA parametrization) and data p4|-[38t. Note that there are several 
other phcnomcnological parametrizations for the deuteron EM form factors in the literature - see Refs. [2^. l25l. . 
Our main observations are: first of all, the one-body diagrams in Figs.3(a)-3(c) correctly reproduce the normalizations 
of charge and magnetic form factors. The contribution of the one-body diagrams in Figs. 3(b) and 3(c), generated by 
gauging the coupling of the deuteron with nucleons is strongly suppressed in all form factors. They become almost 
constant at > 0.5 GeV^. Second, the contribution of the one-body diagram of Fig. 3(a) is dominant in the charge 
and magnetic form factors up to 1 GeV^. Above 1 GeV^ the contribution of the two-body operators Jff^'"^, generating 
the diagrams of Fig.3(d), becomes important for the charge and magnetic form factors. In the case of the quadrupole 
form factor the contribution of the two-body graph [Fig. 3(d)] generated by the operator J^^'^ is dominant for both 
the normalization and the Q^-dependence of this quantity. The operator J^^'^ simulates the contribution of the D- 
wave component of the deuteron wave function. Therefore, inclusion of the two-body operators J'^^'^iq) is sufficient 
to describe the deuteron EM form factors as well as elastic eD scattering. In particular, the two crossing points of 
the charge and magnetic form factors of the deuteron are also successfully reproduced when including the two-body 
operators. Without the two-body operators, one cannot correctly reproduce the two nodes at ~ 0.7 GeV'^(for 
Gc) and 2 GeV^ (for Gm)- 

We also tested different forms of the deuteron correlation function $£> (monopole, dipole. Gauss, etc.). It is 
found that the physical observables are weakly sensitive to its form, whereas they are dominantly controlled by the 
scale parameter of the correlation function. Inclusion of additional terms (e.g. with derivatives) in the Lagrangian, 
describing the bound state structure of the deuteron [see Eq. ([7])], does not lead to a considerable improvement in the 
description of the deuteron observables; it only introduces additional parameters. Because we are interested in keeping 
the number of free parameters to a minimum while still obtaining a good fit to the data we restrict the coupling of 
the deuteron to the constitutent nucleons to the simplest form. 

IV. CONCLUSIONS 

In this work, we applied a relativistic effective Lagrangian approach to study the EM properties of the deuteron 
considering the deuteron as a weakly bound state of a proton and a neutron. We found that in the present approach 
two-body interaction terms are crucial to reproduce the quadrupole moment of the deuteron and the two crossings of 
the charge and magnetic form factors at ^ 0.7 and ~ 2 GeV^, respectively. The effective two-body operators we include 
reflect and model the S/D mixing in the deuteron wave function and pion exchange contributions. With the adjusted 
parameters, listed in Table 1, we obtain a reasonable description of the EM form factors of the deuteron up to ~ 2 
GeV^ using MMD parametrization of the nucleon EM form factors. Note that our result for the deuteron polarization 
tensor T20 using the MMD parametrization deviates from the data and predictions of theoretical approaches [1, [l^, U^l 
in the region ~ 1 — 1.4 GeV^. 

Finally, recent experiments of the electron-proton polarization transfer scattering f4l[ show that two-photon ex- 
change may play a role for determing the charge form factor of the proton [i^j. It is therefore of great interest to 
check if the two-photon exchange mechanism also plays a prominent role on the form factors of the deuteron [i^ as 
well as of the neutron. Work along this line is currently in progress. 
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Table 1. Parameters defining the two-body EM nucleon operators extracted from data using two forms of the 
parametrization of the Dirac and Pauli nucleon form factors: MMD [2lj and Kelly [2^ . 



Ff-(g2) 


Ad (GeV) 


51 (fm^) 


.92 (fm^) 


53 (fm^) 


A (GeV) 


Al (GeV) 


A2 (GeV) 


MMD [21] 


0.25 


0.06 


-1.97 


0.16 


1 


0.70 


0.80 


Kelly [22] 


0.30 


0.23 


-1.48 


0.06 


1 


0.58 


0.58 
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FIG. 4: Form factor | Gc(Q^) I- The solid curve is the result of the TGA parametrization. The double dash-dotted and dashed lines are 
our results with the MMD |2ll | parametrization restricting to one— body and including two— body electromagnetic currents, respectively. 
The double dot-dashed and dotted lines are our results with the Kelly 1231 parametrization restricting to one-body and including two— body 
electromagnetic currents, respectivel y. The data are from [27| (open circle), [28{ (open square), [29l | (open diamond), [30| (Flux), [3]] ] 
(triangle up), .32] (filled circle), and [SSt (filled square). 
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FIG. 6: Form factor | Gm{Q'^) |. The solid curve is the result of the TGA parametrization. The double dash-dotted and dashed lines are 
our results with the MMD ^21] parametrization restricting to one— body and including two— body electromagnetic currents, respectively. 
The double dot-dashed and dotted lines are our results with the Kelly 1231 parametrization restricting to one-body and including two-body 
electromagnetic currents, respectively. The data are quoted from [2^ (circle), [34| (diamond), |35|l (square), |36|| (star). 




FIG. 7: Form factor | Gc(Q^) I- Separate contributions of the diagrams in Fig. 3 with the MMD [2ll| parametrization: diagram in Fig. 3a 
(the dashed line), diagrams in Figs. 3b and 3c (the dotted line), the total contribution of two-body operators generating the diagram Fig. 3d 
(the dot-dashed line). 
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FIG. 8: Form factor | Gq{Q^) \. Separate contributions of the diagrams in Fig. 3 with the MMD [2l| parametrization. 
Notations are the same as in Fig. 7. 




FIG. 9: Form factor | Gm(Q^) |. Separate contributions of the diagrams in Fig. 3 with the MMD [2l| parametrization. 
Notations are the same as in Fig. 7. 
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Q (GeV ) 

FIG. 10: Form factor | Gc(Q^) |- Separate contributions of the two— body operators with the MMD [2l| parametrization: the contribution 
of the jft^'^ operator (the dotted hne), the contribution of the jf^^'^ operator (the dashed line), the contribution of the jff^'^ operator 
(the dot-dashed Une) and the total result (solid line). 




FIG. 11: Form factor Gc{Q'^) on a linear scale. Separate contributions of the two-body operators with the MMD [2l| parametrization. 
Notations are the same as in Fig. 10. 
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FIG. 13: Form factor Gq(Q^) in a linear scale. Separate contributions of the two-body operators with the MMD i21il parametrization. 
Notations are the same as in Fig. 10. 
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FIG. 15: Form factor Gm{Q^) on a linear scale. Separate contributions of the two-body operators with the MMD |2111 parametrization. 
Notations are the same as in Fig. 10. 
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FIG. 16: Form factor A(Q^). The solid curve is the result of the TGA parametrization. The dashed and dotted lines are our results with 
the MMD [2lj and Kelly [22l | parametrizations, respectively. The data are quoted from j29j (circle), [37j (diamond), [36l | (square), [38| 
(star), respectively. 




FIG. 17: Form factor B(Q'^). Notations are same as in Fig. 16 
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FIG. 18: Deuteron polarization tensor T2o(Q^) at de = 70°. The solid curve is the result of the TGA parametrization. The dashed 
and dotted lines are our results with the MMD |21| and Kelly |22|| parametrizations, respectively. The data are quoted from |31| (circle), 
[3I (triangle right), [sl fsguare). J2I (diamond), Mled triangle up), (filled triangle down), (27[ (filled circle), [H (star). The 
theoretical results are taken from Q (dot-dashed line) , [12)| (short-dashed line) , [13|| (double dash-dotted line) , Q (double dotted— dashed 
line). 



